Suppression of Quantum Phase Interference in Molecular Magnets Feg with 

Dipolar-Dipolar Interaction 



Zhi-De Chen^, J.-Q. Liang^, and Shun-Qing Shen^ * 

^Department of Physics and Institute of Modem Condensed Matter Physics, 
Guangzhou University, Guangzhou 510405, China 
^Institute of Theoretical Physics and Department of Physics, 
Shanxi University, Taiyuan, Shanxi 030006, China 
^Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China 

Renormalized tunnel splitting with a finite distribution in the biaxial spin model for molecular 
magnets is obtained by taking into account the dipolar interaction of enviromental spins. Oscillation 
of the resonant tunnel splitting with a transverse magnetic field along the hard axis is smeared by 
the finite distribution which subsequently affects the quantum steps of hysteresis curve evaluated 
in terms of the modified Landau-Zener model of spin fiipping induced by the sweeping field. We 
conclude that the dipolar-dipolar interaction drives decoherence of quantum tunnelling in molcular 
magnets Feg, which explains why the quenching points of tunnel spliting between odd and even 
resonant tunnelling predcited theoretically were not observed experimentally. 



Macroscopic quantum phenomena in magnetic molec- 
ular clusters have been being an attractive field in re- 
cent years[0 |, |, g |, ||, 0, |. Octanuclear iron(III) 
oxo-hydroxo cluster Feg is of special interest because it 
shows not only regular steps in hysteresis curve but also 
oscillation of the tunnel splitting due to the quantum 
phase interference [3]. Oscillation of tunnel splitting of 
the ground state with respect to the external field along 
the hard axis was predicted theoretically by Garg [4] as a 
consequence of the quantum phase interference of tunnel 
paths, and it was subsequently generalized to tunnelling 
at excited states and resonant tunnelling for quantum 
transition between different quantum states with its x 
component of the spin = —10 and 10 — n (along the 
easy axis of Feg) recently [6]. The quenching points be- 
tween even and odd n have a shift 7r/2. However, a se- 
rious problem, why the theoretically predicted shift of 
quenching points of tunnel splitting between odd and 
even resonant tunnelling was not observed in the exper- 
imental hysteresis curves [3,6,9], remains to be solved. 
This is the main motivation of this Letter. Here we use 
the Landau-Zener model [3,8,9,10] to describe the spin 
flipping induced by the sweeping field with a modified 
bare tunnel splitting considering the dipolar interaction 
with environmental spins. There are two basic interac- 
tions to be considered: spin-phonon and spin-spin inter- 
actions. For the molecular magnets Feg in mK tempera- 
ture region, the spin-phonon interaction [8] can be safely 
ignored as the spin-lattice relaxation time is extremely 
long [11]. The interaction between the big spin and the 
environmental spins was considered as the main source 
of decoherence of tunnelling in magnetic macromolecules 
[12] and recently it was shown that the nuclear spin plays 
an important role in magnetic relaxation [13,14]. In this 
Letter, starting from the mean field approximation, the 
dipolar interaction is treated as a local stray field h (see 
the following) with a Gaussian distribution. The tunnel 
splitting in the Landau-Zener transition rate should be 



considered as an average over the local stray field h. In 
doing so we find that the quenching of the tunnelling due 
to quantum interference is suppressed by the local stray 
field, and the steps in the hysteresis curve corresponding 
to odd resonant tunnelling are understood. 

We start with the biaxial spin model for the molecular 
magnets Feg [3-5]. The Hamiltonian is given by [15] 



H = KiSl + K2SI - gugS ■ (B + h), 



(1) 



where Ki > K2 > and B is the external magnetic 
field. The term —gf^^S ■ h is the dipolar-dipolar in- 
teraction between the magnetic molecular cluster and 
the environmental spins, i.e. h = '^jJijSj, where the 
summation runs over the neighboring clusters. Strictly 
speaking, this should be a many-body problem. In this 
Letter, h is treated approximately as a local stray field, 
h — J2j Jij (Sj)- Both experimental [14] and the Monte 
Carlo study[16,17] show that h has a random distribution 
with a distribution width in proportion to (1 — |M|) and 
its mean value proportional to M where M is the total 
magnetization of the system. Here we assume that h has 
a Gaussian distribution with an equal distribution width 
in all directions [18] 
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To simulate the experimental setup [3] , the external mag- 
netic field is taken to be B = {Bx,0, Bz}: a uniform 
field Bz along the hard axis and the sweeping field B^ 
on the easy axis Bj. — nAB ± ct where n is integer, AB 
is the field interval between neighboring resonant tun- 
nelling and c — dB^/dt. In the following calculation, we 
take Ki = 0.310K, K2 = 0.229K, and c = O.lT/sec for 
the molecular magnets Feg [3]. 

Theoretically, quantum tunnelling for a spin system 
without the local stray field can be understood in the 
instanton method [4-6], the Landau-Zener model [8-10], 
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and by diagonalizing the Hamiltonian numerically [3, 
19]. The instanton method can give the tunnel split- 
ting. When the field along the easy axis satisfies the 
resonant condition, + = nAB, the transition rate, 
while the field on the easy axis sweeps over the resonant 
point, is given by the Landau- Zener transition formula 
[8-10], Plz = 1 — cxp (— TrA^/i^n) , where A„ is the tun- 
nel splitting and Vn = 2gijgh(2s—n)c. It should be noted 
that in this way we have assumed tacitly that the tunnel 
splitting for all the spins inside the resonant window are 
the same and thus all the spins tunnel with the same 
transition rate. However, when the local stray field due 
to the dipolar-dipolar interaction is taken into account, 
such a picture should be modified. A random distribu- 
tion of local stray fields like Eq.(2) with a distribution 
width (7 ^ 0.05T typical for the molecular magnets Feg 
will block the resonant tunnelling of either the ground 
state or the low- lying excited states [13]. Nevertheless, 
such a problem can be circumvented by using the sweep- 
ing field along the easy field. When B^ sweeps over the 
resonant point, it will make the spins with different h^'s 
to satisfy the resonant condition, and allows continuous 
relaxation. Since the tunnel splitting is very sensitive to 
the transverse local fields B^ + hz, and hy [4-6], the spins 
tunnel with different tunnel splitting while B^ is sweep- 
ing over the resonant point. Consequently, the spin tran- 
sition rate observed in the experiment should be given 
by 



{Plz) ^ l-exp{-7r(A2 
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where (• • •) represents the average over the distribution of 
the local stray field, i.e., (A^) = / A'^{h)p{h)dh. Accord- 
ingly, the tunnel splitting extracted from the measured 
transition rate should be ^ (A^J but not A„. In other 
words, the starting point to understand the experimen- 
tal observation should be -\/(A^ instead of A„. The two 
quantities are qualitatively different from each other as 
we shall show in the following. 

The instanton method [4-6] is efficient and powerful to 
evaluate the tunnel splitting A„. The Lagrangian for the 
biaxial model Eq.(l) is 



L(n) = —sh{l — cost 
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where |n) is the spin coherent state. With the help of 
the mapping technique, {4>,p = sh cos 9) is regarded as 
a pair of canonical variables. To c;ak;iilate the excited 
state tunnelling or resonant tunnelling, one needs to ap- 
ply the Bohr-Sommcrfeld quantization rule (j) pdcf) = nh 

to define the classical orbits (n is an integer). Then a 
propagator with both imaginary and real time will be 
used to describe the tunnelling between two degenerate 

states, 
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The tunnel splitting is found by integrating over two de- 
generate classical orbits. In molecular magnets Feg the 
local stray field is rather weak, i.e. giJ,^\h\/{K2s) <C 1. 
We calculate the tunnel splitting at the nth resonant tun- 
nelling point with the transverse field Bz + hz, and hy, 
and obtain that, 



A„ « %e-5" {le^'i'^y + e-^^''^ 

+ 2cos[2(s7r- n7r/2 - dnhz - dnBz)]\}^^^ , 
where and 5" is the instanton action 
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V{(1)) = K2S^ sin^ (/) — gfi^nAhs cos (/> 
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En is the energy of the nth excited state, is the turning 
point determined hyV {n—cf)^) = En,Qn is the pre- factor 



Qn ^ A-K/^V"{Q){2Ki+gnBnAh/s), 

q = gnBTT\^'^/2K2{l - A = K2/K^, and 
5Mb r dcj) 



dn = 



2Ki 



^^fi- nAhcosf 



(9) 



(10) 



2K2S 



Using the parameters in Feg, it is found that the contri- 
bution from hy and hz to 5" and thus Qn^~^'' is very 
small under the condition giig\h\/ {K^s) <C 1. The aver- 
age value of A^ is given by 
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where Q„o = Qn{hz = hy = 0), S^q = S'^{hz = hy 
In the absence of the stray field, i.e. a = ho 
above expression reduces to 
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which indicates the oscillation of the tunnel splitting with 
the transverse field and a shift 7r/2 of quenching point be- 
tween the odd and even resonant tunnelling, recovering 
the results in the previous works [3-6]. This is known 
as a result of the quantum interference of the tunnelling 
along two different paths. Qualitative difference between 
^ (A^) and Anihx = hy = 0) can now be seen by com- 
paring Eq.(12) with Eq.(ll). In the case of Bz = and 
integer spin, Eq.(12) predicts that odd n resonant tun- 
nelling quenches due to the quantum interference, while 
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in the presence of the stray field Eq.(ll) gives non-zero 
tunnel splitting 

Vi^, ^ ^Qnoe-^-y^V^^T^=^ (13) 

for ho = 0. The quenching due to the quantum inter- 
ference is suppressed by the local stray field. In another 
word the quantum tunnelling for odd n is decoherenced 
because of the dipolar interaction with the environmental 
spins. The tunnel splitting of all six resonant tunnelling 
for the molecular magnets Fcg with and without the lo- 
cal stray field are shown in Table I. We see that y^A^J 
for an odd n increases from zero while the random field 
become stronger. The random field also increases the 
tunnel splitting of even resonant tunnelling. It increases 
about 2.7 times as a becomes as large as 0.08T, which 
resolves the puzzling that the experimental observation is 
about 3.0 times larger than the numerical result for the 
tunnel splitting [3]. A detailed evolution of the tunnel 
splitting with the distribution width around the topolog- 
ical quenching points is shown in Fig.l. As the width 
of the distribution is proportional to (1 — |M|), the cal- 
culated results for different M arc shown in Fig.l, which 
are in good agreement with the experimental observation 
(see Fig.lO in Ref.[20]). One can see from Eq.(ll) that 
the main effect of ha is to provide an initial phase and 
thus shifts the oscillation. For Fcs, ho — cr/4 [16], and the 
effect of modification for nonzero ho is almost omissible. 

Table I: Trmnel splitting -y/ (A^ ) (the unit is Kelvin) 
for Fes in the case of = ho = 0. 



FIG. 1: Illustration of ^/{^^ (n = 0, 1) around topological 
quenching points due to quantum interference with different 
distribution width a's. 
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The oscillation of the tunnel splitting for ,s = 10 with 
various distribution width a's is shown in Fig. 2. From 
Fig. 2, it is shown that the oscillation of the tunnel split- 
ting due to quantum interference is suppressed by the 
local stray field h. For a distribution width a = 0.05T 
which is estimated for Fes [9,13], the oscillation of tun- 
nel splitting with respect to the field along the hard axis 
is still visible, while the oscillation is suppressed com- 
pletely for the width as large as 0.08T. In fact, when the 
distribution width approaches the half oscillation period, 
the oscillation due to quantum interference disappears 
and the classical behavior, i.e. tunnel splitting increases 
monotonously with Bz, is resumed. The above analysis 
leads to a decoherence mechanism for quantum interfer- 



FIG. 2: The oscillation of y(A§) with different distribution 
width a's for s=10. Prom top to bottom: a = 0.08T, 0.05T, 
0.02T, and O.OT. 

ence due to the dipolar-dipolar interactions between the 
spins without dissipation [12]. 

The magnetization jiunp from the spin flipping at the 
resonant tunnelling can be calculated from the modified 
Landau-Zener transition rate given in Eq.(3). In princi- 
ple the time evolution of the spin system in Eq.(l) can 
be obtained by solving the time-dependent Schrodinger 
equation ih-^ |<I>) = |$), which contains a set of (2s-|-l) 
coupled differential equations for the model in Eq. (1). 
It was shown [19] that the coupled differential equations 
can be reduced to that of an effective two-level system 
with the effective Hamiltonian. Here we have 

f -{W-n)gi,sCt ^/(M/2\ 
Hesit) = , (14) 

V v^(A2)/2 lOffMBCt / 

near the resonant condition and the time-dependent state 

isgivenby |$cff) = a_io(^) |-10)+ aio-„ (0 [10 - n). The 
tunnelling splitting in Eq.(14) is ^ (A^) instead of A„ as 
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FIG. 3: Hysteresis curves with different distribution widtfi 



we discussed. Correspondingly, the magnetization jump 
from the nth resonant tunneUing is obtained as 

AM„ = ($eff| |$off) |t=+oo - ($cff| S.^ |$cff)U_o„ . 

(15) 

Numerical results are shown in Fig. 3. It is worth em- 
phasizing that the resulting jump is modified as a rather 
smooth one due to the local stray field. The hysteresis 
curves in Fig. 3 are drawn with the initial condition of 
Sx = —10, i.e.,a_io(i = — oo) = 1. As it is shown in 
Fig. 3 the steps in hysteresis curve are smeared grad- 
ually with increasing the distribution width of the local 
stray field and the observed curve in experiment [3,9] can 
be recovered from the present theory. 

In this Letter, the local stray field is treated as a 
"frozen" one inside the resonant window. Strictly speak- 
ing, both the width and the mean value of the distri- 
bution of the field should vary with the time-dependent 
magnetization during the resonant tunnelling. However, 
it should be noted that a "frozen" distribution is based on 
validity of the Landau-Zener model. If the spins "feel" 
the change of the local field due to spin flipping, the 
spin transition rate is no longer the one in Eq.(3). In 
that case, one should consider the non-linear Landau- 
Zener tunnelling [21] and the "hole-digging" mechanism 
[13,14]. This indicates that our result is valid when the 
field sweeping rate is not too small such that the evolu- 
tion of the local field is relatively slower than the sweep- 
ing field. Namely, the overlap time of two levels in res- 
onance Ti ~ A„/(2/Lt^5'c) should be less than the char- 
acteristic relaxation time T2 due to the dipolar-dipolar 
interaction. This means that the field sweeping rate 
c > Co ~ An/{2fj.gScT2)- In Fcg [14,16] the ground state 
tunnelling Aq ^ 10^ ''K, T2 ~ lO^^sec, and Cq is esti- 
mated to be 10~^T/sec., which is in good agreement with 
the experimental condition [3,20]. On the other hand, a 
finite distribution of tunnel splitting due to the local stray 
field has a deeper impact on the magnetic relaxation. If 
all the spins tunnel with the same tunnelling rate, the 



magnetic relaxation should obey the exponential law, i.e. 
e~'"* where F = 2Plzc/A where A is amplitude of the ac 
field used in the experiment [20]. Instead, in the present 
picture, there is a finite distribution of tunnel splitting 
p{An) which will lead to a finite distribution of the re- 
laxation rate p{T) characteristic of the complex system 
like spin glass [22]. Consequently the resulting relaxation 
will obviously deviate from the simple exponential law as 
observed in the experiment [20] . Further analysis will be 
provided elsewhere. 

We have studied the effect of dipolar interaction be- 
tween giant spins in the molecular magnets Feg in the 
mean field approximation which leads to a Zimman term 
of the spin in the local stray field. Our main observation 
is that the topological quench due to the quantum phase 
interference of tunnel paths is suppressed by the finite 
distribution of the local stray field, and the steps in the 
hysteresis curve corresponding to odd resonant tunnelling 
are explained theoretically. Thus we conclude that the 
dipolar-dipolar interaction leads to the decoherence of 
quantum tunnelling in Fcg. Finally it is worth pointing 
out that the mechanism of decoherence may not be just 
limited in Feg, but can be generalized to other molecu- 
lar magnets such as Mni2 since the local stray field due 
to the dipolar-dipolar and hyperfine interactions always 
exists. 
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